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Abstract 



The effect of projectile shape on cross sections and momentum distribu- 
tions of fragments from heavy ion reactions is studied. We propose a new 
approach that implements the underlying symmetries of each isotope with a 
few parameters directly in the density. Various densities and their nuclear 
structure are then analyzed in the reactions of ^^C and ^^Li, ^^Be, and ^^C 
on a carbon target. 
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I. INTRODUCTION 



In a previous paper [|l| in reactions of ^^C and ^^Li projectiles on a ^^C target, it has been 
shown that the total reaction cross section, inclusive fragment cross sections and fragment 
momentum distributions are sensitive to the assumed proton and neutron radial density 
distributions of the projectile. These results stimulated us to further investigate effects of 
nuclear shape of neutron rich isotopes in the heavy ions reactions. A special attention in 
our investigations will be given to the neutron halo nuclei for the following reason. Recent 
developments of radioactive nuclear beams enabled one to study a detailed structure of 
nuclei far from the stability line [H. The neutron halos have been observed in nuclei near the 
neutron drip line by reaction measurements with intermediate and high energy radioactive 
nuclear beams. These nuclei are of a special interest in relation to shell structures near the 
drip line and to new excitation modes associated with the excess neutron on the nuclear 
surface. 

In this paper we use the method of hyperspherical functions (HSF) to construct the 
densities of the light drip line nuclei. Such nuclei typically have small binding energies 
and extended radial densities, and for them a suitable treatment of the tails of nuclear 
radial wave functions is essential. Physically such wave functions should have exponential 
tails. A main advantage of the HSF is that it provides a realistic radial wave functions 
at large distances. In addition, the spurious center-of-mass motion is removed from the 
onset and the symmetries are properly taken into account. In [Q, an attempt was made 
to provide a unified description of ^''''^'^'i^Li rather than just a single isotope. Instead of 
trying to parameterize the effective interaction for each isotope, a simple parametrization 
for all the isotopes is used. Furthermore, there is no inert core and all the nucleons are 
properly antisymmetrized. Finally, because Jacobian coordinates are used, no problems are 
encountered with the treatment of the center of mass. The binding energy of a nucleus 
serves as a natural scale in this approach and is related to all other observables. 

In order to account for the clusterization effect more precisely in our calculations, we ex- 
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pand the HSF in the translationally invariant many body harmonic oscillator (HO) functions 
1^ and implement the nuclear structure straight in the densities. We use a few one-body 
oscillator parameters in order to better reproduce the experimentally observed mean-square 
radii (rms) and density distributions. These densities are then used in Monte Carlo sim- 
ulations of nucleon-nucleon collisions, in which the total and inclusive cross sections and 
transverse momentum distributions of fragments are calculated. Thoroughly tested, this 
approach gives a clear picture of nuclear shape. 

The paper is organized as follows. In Section II, the HSF method is described. Next 
we elaborate on the details of the expansion with respect to the radial wave functions of a 
(3A-3)-dimensional harmonic oscillator. In Section HI, we give a brief summary of the light 
isotopes and construct their Young diagrams. Knowing the underlying symmetries, we first 
build the "normal" ^^C density and compare it with other approximations. We calculate 
the elastic form-factor and run a ^^C + ^^C simulation to see how it works in a standard 
case. Armed with the results, we proceed to the neutron rich ^^Li and also ^^Be, ^^C for 
comparison. In order to make the paper more coherent, we include all the densities in Section 
HI and analyze the reactions in Section IV. The results are in reasonable agreement with 
experiment and other theoretical approaches. Our conclusions are summarized in section V. 



II. METHOD OF THE HYPERSPHERICAL FUNCTIONS 

The HSF method involves a collective variable (hyperradius p) which is related to the 
mean-square radius of the nucleus p'^ = A < >, i.e., to the mean nuclear density. The 
excitations of this degree of freedom correspond to the monopole oscillation of the nucleus 
as a whole and thus the density can be treated as dynamic variable. 

The wave function \I' of a nucleus with A nucleons is translationally invariant and depends 
on the xi,X2, ■ ■ ■ , xa^i Jacobi coordinates, defined as 

xi = W^(ri -ra) 
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X2 



-(ri + rs) -ra 



(1) 



lA-1 



XA-1 



A-1 



^ i=l 



rA 



Here Tj is the coordinate of the i— th nucleon. Each coordinate is a distance between the 
{i + l)-th nucleon and and the center of mass of the groups of nucleons 1,2, . . . , i. The 
hyperspherical coordinates (hyperradius p and {3A — 4) hyperspherical angles) are chosen 



as 



xi — p sin 9n-i ... sin 02 sin 9i 
X2 = p sin 9n-i ... sin 6*2 cos 9i 



Xn-l = psinOn-l COs9n-2 
Xn = PC0S6'„_1 



(2) 



< p < oo 



1=1 



0< 9i<2n, 0<9k<n k>l 



for n — 3{A — l)-dimensional space of the Jacobi coordinates. For this set of coordinates, 
the volume element dV reads 

dV = dxi dx2 ■ ■ ■ dxn = p"^^ dpdVL (3) 

where the solid angle element dO, is 

dn = sin"-2 9n-i sin"-^ 9n-2 ... sin 92 d9n-id9n-2 ■■■d9i (4) 

The Laplacian is given by 

"^"^ = sin-'^„_,a^„_, + il^^""- 
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The hyperspherical functions, or i^— harmonics, are the eigenfunctions of the angular 
part of the Laplacian 

^n^YK^ie^) = -K{K + n- 2)F,,^(^,). (6) 

The value of K is the analog of the angular momentum at n = 3 and is called the global 
momentum. The subscript 7 denotes all the quantum numbers necessary to enumerate the 
various degenerate states of Eq. (||). For 7 it is expedient to use the Young diagram [/] 
and the Yamanouchi symbol (r) , characterizing the /T— harmonics properties relative to the 
A— nucleon permutations, and LST to designate the orbital momentum, spin and isospin of 
this state. 

The wave function of a nucleus is then expressed in the form of an expansion in the 
i^— harmonic polynomials 

K-y 

POO 

/ xl,{p)dp = l l = [f]rLST (7) 

The Hamiltonian reads 

2m p3A-4 dp \ dp) 2m ^ ^ 

Here m is the mass of a nucleon. The Schroedinger equation for the radial function XK-yip) 
can be written as 

(;^-^^^^^%^-|?[^ + <>)]|xx>) = |f E <^'(P)XW(P) (9) 
{ap p n ) n ^,y_^^^ 

where Lk = K+{3A- 6)/2, and VT^^'^ (p) are the matrix elements of the nucleon-nucleon 
interaction 

v = j:yinj) vin,) = fiu,)w^, (10) 

i<j 

which may be expressed through the two-body fractional parentage coefficients in the form 
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Wk^^\p) = < AK[f]rLSTMLMsMT\V\AK'[fy L' S'T' Ml'Ms'Mt' > 

= ^^^Y^ E < AK[f]rLST\iA - 2)K,[f]r,L2S2T2AiL" K"); LoSoTo > 
X < SoTo\W^r\SoTo > R^ff^P) 

X < (A - 2)K2[f]r2L2S2T2A{L"K"y, LoSoTo\AK'W L' S'T' > (11) 

Here 

Rk-lM = J dOiSm^^-^ei COS^ eiUKK"Lo^K'K"Lo 

xPx^i"^^/'"'''^''^'/'(cos2^i) (12) 

< S'oTolP^o-T-IS'oTo > is the spin-isospin part of the interaction matrix element. The collective 
potential (p) is given by integrals containing two-body matrix elements together with 

the i^'— harmonics fractional parentage coefficients. In the previous investigations [0 it 
was shown that with the HSF one may use the same formulas for the fractional parentage 
coefficients which were obtained earlier with the translationally invariant shell model. In this 
case the Talmi-Moshinsky coefficients must be replaced by the Raynal-Revai coefficients, 6j- 
symbols must be added for the over-binding of the global momentum K, and additional phase 
multipliers are inserted. Their use has led to a significant improvement of the computations 
P). The calculation is simplified if only the first few terms in the expansion i.e., K = 
Kmin and K = Kmm + ^ are taken into account. Usually one adopts the Kmin approximation 
in which all values of K greater than Kmin = A — 4 {A < 16) are neglected. The success of 
this approximation lies in the fact that the centrifugal barrier reduces the contributions of 
configurations with K greater than Kmin in the equations determining the hyperspherical 
wave functions. 

A. Expansion 

The system has an analytical solution for the harmonic oscillator potential V = 
mc(j^p^/2 0. The energy eigenvalues are given by 
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Tl 

S, = {2u + K+ -) huj (13) 



with the eigenfunctions 



2^^' .e-^^^x^Lf^^("-^\a:2) (14) 



^r(z/ + ir+f) 

where we have introduced a dimensionless parameter x = p/r^ with Tq = h/muo and the 
associated Laguerre polynomials are defined as 

Let us briefiy examine critical regimes for the Schroedinger equation (P). As p — > 0, it is 
reduced to 

^Xip) = Xip) (16) 

with the trivial eigenfunction 

X{P) - P^^'+' (17) 

On the other hand, as p — > oo, the Eq. (||) reads 

d'^ . . 2mE . . 

5^X(P) = ^Xif) (18) 
and the eigenfunction has the following asymptotical behavior 



X(p) ~ exp I -y pj (19) 

as was mentioned in the Introduction. We will return to this important point in next the 
section. 

Now we can express radial functions xkip) terms of the 3(^4 — l)-dimensional harmonic 
oscillator, as shown on Figure 1 for the ^^C nucleus, 

xI{p) = T.cIkR.k{p) (20) 

V 

and the expansion coefficients are given by 
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CIk = J XI{P)RMP) P^^'^-'^ dp (21) 
Co ~ 0.984, Ci ~ -0.061, 
0.158, Cg- 0.055. 

with the oscillator frequency hu determined so that the maximum of the lowest oscillator 
function Rok{p) coincides with the maximum of the lowest hyperspherical function Xic(p)- 
We will deal only with the ground state wave function and use C^Kmi„ ~ short. 

It is striking that the overlap integral of the ground state radial function Xxip) with 
the lowest oscillator function Rqk{p) is 98%. At the same time, the contribution to Xxip) 
of the two-quantum oscillator excitation is exceptionally small. A very small value of the 
coefficient Ci is not fortuitous but conforms to the theorem |^, which states that if the 
oscillator frequency hu is chosen to make the overlap integral of Rqk{.p) and Xxip) maximal, 
then the overlap integral of Xxip) with Rik{p) is strictly zero. 

As we have already noted, in choosing hu we did not achieve a maximum value of the 
Xxip)^ Rok{p) overlap, but assured only a coincidence of the maxima of these functions. 
However, if such a condition is satisfied, the overlap integral of the functions is fairly near the 
maximum, and the coefficient Ci is very small, although nonvanishing. This circumstance 
casts light on the origin of close agreement between the results of the translationally invariant 
shell model and the HSF method for the ground state. Since the function Rik{p) cannot 
be mixed with Rqk{p), the shell function in the HSF method is improved by admixture of 
states with oscillator energy Ahu and higher. Because of the large difference between the 
energies of these states and the ground state, this admixture is small. 

Then one can derive the binding energy as a function of the interaction matrix elements 

1. 5A-11. 



2 ' 2 

Cu+i L \ , 5A-11, , \ , , 5A-13. 



+ 2^^ 



+ [ KK{p)w^i:{p)R,K{p)p'^-'dp (22) 
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The rms of a nucleus is given by 

.c; 4 - 1 1 I .c; 4 - 1 1 

(23) 



h 1 

< >. = ^ 

muj A 



III. SYMMETRIES OF NEUTRON RICH ISOTOPES 



Using the Young- Yamanouchi classification scheme fly] , which stems from the orthogonal 
representation of a given permutation group, and with the knowledge of the total spin and 
isospin of the nucleus, we determine the symmetry with a corresponding Young diagram 
[/] of each isotope. In Table I results of the classification of ^"^'^Li, ^~^^Be, ^"^'^B (for 
completeness), and ^-^^C nuclei are collected. The Lithium sequence, illustrated in Figure 
2, reflects the strong excitation in the isospin. The neutron super-rich ^^Li obtains an extra 
pair of halos in the same manner as ^^Li results from ^Li by adding two halo neutrons. In 
the brackets we have put two possible versions of a transitory ^°Li state. 

Note that while for the most part these sequences are rather straightforward, obtained 
by a simple adding of a neutron to the p-shell, some isotopes like ^^Be do stand out. There 
has been a great deal of interest in the effect of clusterization, and our scheme seems to 



agree with a recent study of the multi-cluster structure that a simple description in 
terms of one Young diagram may not be adequate for some isotopes. 

Having thus established the symmetry requirements on the neutron content of different 
isotopes, we proceed to construct the density distributions for the light nuclei. Originally, 
this classification was applied to the construction of the HSF wave functions In Table 
II we reproduce their results for the binding energy and rms for ^~^^Li and ^^C isotopes in 
comparison with the experimental data. It can be seen that the variation in the binding 
energy as a function of the mass number and rms for the Lithium isotopes is qualitatively 
reproduced. A comparison of ^^Li and ^^C shows a very strong nuclear structure effect on 
the binding energy and rms. 

If one would like to look at the more detailed description of the nuclear structure which 
comes with the wave function and exhibits the nuclear shape through the nuclear densities, 
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one may use the information that comes from measurements of the total reaction cross 
sections and the inclusive cross sections for fragments in which the projectile has lost one 
neutron or proton. In Ref. |T| it has been shown that both the cross sections and the 
fragment momentum distributions lead to the conclusion that the HSF density distributions 
in the -ft'mm approximation extend too far for protons and not far enough for neutrons. The 
inclusion of several more terms in the expansion (J^) will give an extension of the asymptotic 
density to larger radii. This will give better agreement between the HSF method and the 
experimental data. For this reason we propose another approach in which the clusterization 
effect is taken into account directly in the density. 



A. Density 



From the formulae in section II one derives the density of nuclear matter through 



K 



J 



/7r' 



with 



n 



and the radial component is expressed as 



(24) 



(25) 



J" 



(26) 



The radial density distribution reads {A < 16) 

16 r(i(5A-ll)) .oo (p2_^2)|(5A-16) 

= /-^ n . 7TTY / :^73T^ X^{p)xAp) dp 



v^r(i(5A-14) 



A-4r(i(5A-ll) .oo ,2(^2_^2)^(5A-18) 

+ - — 7= 7T— rv / ^T^. XiKPjXAP) dp 



3 r(i(5A-16 



P 



(27) 



with the s-shell (first line) and p-shell densities explicitly shown. We also treat the halo 
neutrons Nh separately from the p-shell nucleons. The rms of each shell is calculated as 
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2 ^ / n^^^'-'-(r)r^dr 



<ru> shell j^sheii^^y^^ (28) 
with the density normahzation 

47r r n^'^'''^ {rydr = {4, A - 4 - Nh, Nh ) (29) 
Jo 

Before calculating the densities for the neutron rich isotopes, we would like to check on 
a standard nucleus such as ^^C, for which there exist a few other density approximations. 

• In Ref. the exact HSF density for ^He, ^Li, ^^C, and ^^O was approximated by 

nir) = E a. 5 e-^Vr? (30) 



i=0 

motivated by the analogy with the harmonic oscillator. Using different oscillator pa- 
rameters Cj, Ti, one could easily fit various experimental observables, such as the rms 
and so on. For the ground state of ^^C, it was enough to have the first two terms 

n(r) = ase-''"/'' + flp ^ e""'/'^^ (31) 

with Tg = 2.0 fm, = 1.7fm and = 0.718348, = 0.957798 fixed to the normal- 
ization of 4 s-shell and 8 p-shell nucleons. 

Harvey used a phenomenological approach |13| to obtain the densities for ^^C and ^^Li. 



Since the Monte Carlo simulations require a separate density distribution for nucleons 
with different binding energies, a harmonic oscillator model was used for projectiles 
and a standard Fermi shape for the ^^C target. For example, for ^^Li the s— shell 
parameter was chosen to be 1.619 fm for protons and neutrons, and for the p-shell 
2.0 fm for protons and 2.1 fm for neutrons. It was also assumed that there was one 
proton and four neutrons in the p3/2-shell. The outer two neutrons were assumed to 
have an exponential density distribution beyond a radius of 2.5 fm 

-2r/L 

n(r) = 0.3323^;^ (32) 
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with L=3.7 fm. As we mentioned above, the HSF functions have an exponential 
asymptotics ( ]T9D and thus become very attractive from a phenomenological point of 
view. The normahzation, 0.3323, came from the requirement that the volume integral 
contain just two neutrons. The values of the harmonic oscillator parameters and of L 
were chosen to fit the experimental reaction cross section and the ^Li, ^He inclusive 
cross sections for the reaction ^^Li+C at 790 MeV/A. 

• Returning to the HSF method, one does not have as much freedom in choosing various 
parameters. It is very important that in the HSF one obtained a self-adjusted system 
in which the parameters of the interaction determine the size of the system. How- 
ever, working in the -ft'mm approximation, one needs to improve the neutron density 
asymptotics. Having derived the explicit formula (^) for the density distribution and 
knowing the symmetry requirements on the neutron content of different isotopes, we 
are now in the position to do it directly in the density. 

Since the wave function characterizes the nucleus as a whole, the binding energy Eqk, 
cf. Eq. (0), sets its scale = h/mu. It also serves as a natural parameter for the s-shell, 
the core of a nucleus. The p-shell basically carries all the information about a given isotope, 
and its scale will be chosen so that experimental rms is reproduced. In accordance with 
our group-theoretical framework, another parameter will be used for the halo neutrons, still 
fitting to the experimental rms. 

For example, the p-shell density then reads 

|(5A-18) 



"^^^^ = 3^r(i(5A-16)) Jr X.iP.)xAPA) dp, (33) 



with scaling 



Pp^ — PA^ — (34) 



Table HI shows the parameters for these densities, the rms for each shell and of the 
nucleus as a whole. They reproduce the experimental radii of Tanihata, et al. In 
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Figure 3 we plot the densities for ^^C and ^^Li. The ^^Li one is in reasonable agreement 
with the shaded region of the experimental density p5| . 



B. Charged form- factors 

Having obtained the total densities, one can quickly probe its charge distribution by 
looking at the elastic e~ scattering off that nucleus, or its charged form-factor. For a nucleus 
with Z protons it may be written in the Born approximation as 

FcK{q) = fM fcM{q)\FM) (35) 
with the elastic point-proton form-factor Fzifl) defined as 

(g) = 47r [rf^^^r'^dT (36) 

JO qr 

where Uzir^ is the proton part of the total density. The correction for the center of mass 
motion jcuifl) is not required in our approach since it is factorized from the very beginning 
(cf. Eq. 1^). /p(g) is the correction for the finite proton size, which we take to be 

/p(?) = [l + ?V%']-' (37) 

with go = 4.33 fm~^ corresponding to r^ms = 0.8 fm for the proton. 

In Figure 4 we plot the charged form-factors for ^Li and ^^C with the data points taken 
from ||l6|,|l7l, where for the momentum axis a standard ge// is used 



9eff = qn + i^—^] (38) 

rms-^0 

Note, that the elastic rms radii for these nuclei are different from the reaction ones 0, 
and we have slightly modified our parameters so that their densities reproduce the observed 
rrm,s(^Li) = 2.56 fm and r^msi^'^C) = 2.48 fm. Our results match well the minima of the 
experimental form-factors and qualitatively follow their shape. They also compare favorably 
with the full HSF calculation 

Having determined all the parameters in the density, we are ready to analyze the reac- 
tions. 
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IV. MONTE CARLO AND NUCLEAR SHAPE 



The Monte Carlo method allows the calculation of cross sections and momentum distri- 
butions of primary fragments arising from nucleon-nucleon collisions between a target and 
a projectile at energies above 200 MeV/nucleon. The primary fragments may subsequently 
decay to the final observed nuclei. 

This method has been shown |TP| to give results in good agreement with experiment 



for the total reaction cross section, fragment inclusive cross sections, fragment exclusive 
(coincidence channel) yields, fragment momenta, both inclusive and in different coincidence 
channels. It is also possible to calculate the momentum distributions of nucleons that are 
either knocked out of the projectile or come from the subsequent decay. Thus, unlike many 
theoretical techniques, the Monte Carlo decay method permits the comparison of calculated 
and experimental values for essentially all the measured quantities and allows to test radial 
density distributions derived from various theoretical models. The purpose of the present 
work is therefore to compare experimental values with those obtained from the HSF density 
distributions with the aim to get the maximum of information about nuclear shape. 



A. ''C + '^C 

When one or two neutrons are lost in a collision between the projectile with a diffuse 
neutron rich surface and a target nucleus, the transverse momentum distribution of the frag- 
ments thus produced can be fit with a double Gaussian curve in which a narrow component 
is superimposed on a much broader one [|l3l. This phenomenon is also observed for frag- 



ments from the normal projectile ^^C. The shape of the fragment momentum distribution 
is sensitive to the assumed projectile surface diffusivity. As the surface diffusivity increases, 
the width of the broad component decreases. For the sharp cutoff shape, the distribution 
has no narrow component. 

The study of the nuclear structure effect in ^^C + ^^C scattering, using the optical 
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model and the calculation of the total reaction cross sections in the Glauber model 



|2T| convincingly show that the total reaction cross sections are sensitive to the form of the 
density. With that in mind, we have conducted the Monte Carlo simulations with three 
different densities: minimal n{RoK), obtained with Rok instead of the full function x of Eq. 
(PPP, n{x) and the approximation of Eq. (0). The experimental data and our results are 
collected in Table IV. As expected, the total cross sections did not vary much, since we had 

the S3jIX1G Tins 

for each density distribution. 

The tail effect is deduced from the momentum distribution of ^^C fragments. In a fit 
with the double Gaussian curve, the narrow component reflects the behavior of the tail of 
the density and thus we see that having all four functions of the Eq. (pO]) decomposition 
improves the resulting width of both narrow and broad components as compared to the 
minimal ti^Rqk) density. 

In Figure 5 we fit the transverse momentum distribution of the ^^C fragments from our 
n{x) density with one and two Gaussian curves. One Gaussian fit is quite reasonable and 
gives a very close (Tb = 103 to the experimental cr^ = 105. Overall, our numbers match well 
with the experimental data and phenomenological results of Ref. [ffB . 



Having thus convinced ourselves that this approach works well for the standard nucleus, 
we turn to the reactions with the neutron rich isotopes as projectiles on the carbon target. 

B. "Li + 

In analyzing this reaction, we emphasize the halo effect. The momentum distribution 
of ^Li fragments from ^^Li + C at 790 MeV/c has been experimentally studied p2|. It 



has broad and narrow Gaussian components with widths of 80 and 21 MeV/c and a ratio of 
areas B/N of 1.5. Since the shape of the momentum distribution is sensitive to the projectile 
surface diffusivity, we have used two different HSF densities with a very loose pair of halo 
neutrons ni^^ljih) and another one ni^^ljip) by placing them relatively close to the p-shell. 
The results for ni^^'Lih) are shown in Figure 6 and the rest is summarized in Table V. While 
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the width of the broad component is well reproduced, the narrow component is too wide, 
especially for r;,(^^Lip). This result is consistent with the absence of a long neutron 'tail' 
at large radii and low local neutron densities which give a small value for the local Fermi 
momentum. A short tail also manifests itself in the somewhat low total cross sections, which 
drops even more for n(^^Lip). This indicates that the halo pair gives a major contribution 
to the neutron tail and has an unusually large rms radius. We would like to point out that 
our result r"^ >o ~ 5.23 fm for an outer pair of halo neutrons for n(^^Lih) is consistent 
with experimental 4.8 ± 0.8 fm [|l^ and 5.llo;gfm . 



At last, we would like to analyze the structure effect by comparing ^^Li, ^^Be, and 
^^C, which have the same number of nucleons but different Young diagrams. From ^^C 
with [/] = [443] and isospin T = 1/2, one sees the increase in the isospin in ^^Be with 
[/] = [4331], T = 3/2 to i^Li with [/] = [4322], T = 5/2. Their cross sections, see Table 
VI, vividly show this change in their nuclear structure. Starting with auc = 880 mb, the 
cross section jumps to cr^ = 946 mb for ^^Be, indicating some sort of a halo neutron tail and 
culminates in au^i = 994 mb with two halo neutrons. It unambiguously verifies that our 
approach correctly implements the underlying symmetries of each isotope and works well 
with a judicious choice of parameters directly in the density level, bypassing the parameters 
of the nucleon-nucleon interaction. 



V. CONCLUSION 

The hyperspherical functions method provides a convenient basis for the description 
of the neutron rich isotopes. The underlying symmetry of each isotope exhibits a simple 
structure and reproduces the clusterization effect. Studies of the reaction cross section and 
the inclusive cross sections for fragments in which the projectile has lost one neutron or 
one proton give a powerful method for investigating the neutron, proton and total densities 
in the surface of the projectile. Here we have proposed a new approach which takes into 
account the clusterization effect directly to the densities of the drip line nuclei. 
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TABLES 



TABLE I. Classification of the isotopes 







T 


[f] 


6Li 


1+ 





[42] 


^Li 


3/2- 


1/2 


[43] [421] 


8Li 


2+ 


1 


[431] [422] 


9Li 


3/2- 


3/2 


[432] [4221] 




(1+) 


2 


[4222] 




(2+) 


2 


[4321] 


iiLi 


3/2- 


5/2 


[4322] 


i3Li 


3/2- 


7/2 


[43222] 


6 Be 


0+ 


1 


[42] 


^Be 


3/2- 


1/2 


[43] [421] 


8Be 


0+ 





[44] [422] 


9Be 


3/2- 


1/2 


[441] [432] [4221] 


lOBe 


0+ 


1 


[442] [4321] [4222] 


iiBe 


1/2+ 


3/2 


[4331] 


i^Be 


0+ 


2 


[4431] [4422] 


7B 


3/2- 


3/2 


[421] 


8B 


2+ 


1 


[431] [422] 


9B 


3/2- 


1/2 


[441] [432] [4221] 




3+ 





[442] [4321] [4222] 


iiB 


3/2- 


1/2 


[443] [4421] [4322] 


12b 


1+ 


1 


[4431] [4422] 


13b 


3/2- 


3/2 


[4432] 


8c 


0+ 


2 


[422] 




3/2- 


3/2 


[432] [4221] 


10c 


0+ 


1 


[442] [4321] [4222] 


lie 


3/2- 


1/2 


[443] [4421] [4322] 


12c 


0+ 





[444] [4422] 


13c 


1/2- 


1/2 


[4441] [4432] 


14c 


0+ 


1 


[4442] 
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TABLE II. Calculated binding energies and the rms radii. 



A n 


huj [Mev] 


- Eo [MeV] 


- Ee^p [MeV] 


r p 1 
Trms [imj 


1.2 A^'-^ [fm] 




19.47 


30.1 


30.5 


2.41 


2.18 


7-r ■ 


19.71 


38.0 


37.7 


2.45 


2.30 


8T,i 






3Q 7 


9 57 


9 40 


9Li 


16.21 


42.4 


43.8 


2.63 


2.50 


iiLi 


13.68 


48.1 


44.1 


2.76 


2.67 


HQ 


16.01 


72.3 


70.4 


2.61 


2.67 



TABLE III. Parameters and rms radii for each shell, in fm. stands for the parameters of the 
outer/halo neutrons. ^^Lip has four of them close to the p-shell; for ^^Li/j two neutrons are on the 
p-shell, and another two are loosely bound halos. 







r-p 


ro 


V< r'^ >s 


y< r2 >p 


V< r'^ >o 


V< r'^ >tot 


6Li 


1.46 


1.52 




1.84 


2.50 




2.09 


"Lip 


1.74 


1.93 


2.77 


2.18 


2.95 


4.48 


3.14 


"Li, 


1.74 


1.93 


1.93 


2.18 


2.95 


2.95 


3.14 


halo 






3.28 






5.23 




"Be 


1.64 


1.75 


2.63 


2.05 


2.81 


4.26 


2.86 


lie 


1.60 


1.71 




2.01 


2.76 




2.52 


12c 


1.61 


1.77 




2.02 


2.86 




2.61 


^'C Eq. (|1D 


2.0 


1.70 




2.45 


2.68 




2.61 



20 



TABLE IV. Experimental and calculated total reaction cross section (tr and inclusive cross 
section au^ for the reaction ^^C + C, 1000 MeV/A. ub and are the widths of the broad and 
narrow momentum components. Cross sections are in mb. 







a N 






Experiment 


105 




930 


47 


Ref. 


96 


36 


929 


61 


n{Rm) 


87 


53 


943 


65 


n{x) 


96 


40 


937 


63 


n(EqJl) 


113 


45 


929 


48 



TABLE V. Comparison of transverse momentum widths for ^Li fragments, experimental and 
calculated total reaction cross section and inclusive cross sections for the reaction ^^Li + C, 790 
MeV/A. aB and are the widths of the broad and narrow momentum components. B/N is the 
ratio of their areas. Cross sections are in mb. 







(J AT 


B/N 










Experiment 


80 


21 


1.5 


1042 


213 


62 


26 


Ref. 


80 


21 


1.22 


1036 


180 


54 


30 


n(iiLi,) 


80 


28 


1.65 


994 


147 


48 


30 


n(iiLip) 


82 


34 


1.75 


978 


147 


37 
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TABLE VI. Comparison of Li, Be and C projectiles on a standard C target. All 
cross-sections are in mb. 









Ref. 1 


241 


"Li 
"Be 
"C 


994 
946 
884 


1042 
942 


800.6 
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Figure 2: Young Diagrams for Lithium isotopes. 
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Figure 3: (a) The HSF density distribution for ^^C (solid line) as compared to the Eq. (31) 
approximation (dashed) and phenomenological one (dot-dashed) from Ref. [13]. (b) The HSF 
density for ^^Li. 




Figure 4: The charged form-factors for (a) ^^C and (b) ^Li. Experimental points are taken from 
Refs. [16] and [17] respectively. 




transverse momentum C (MeV) 



Figure 5: (a) Transverse momentum distribution of ^^C fragments from the reaction ^^C + C, 1 
GeV/nucleon. Filled squares are the calculated points, with a bin size 20 MeV/c. Dashed lines 
are least squares fit to the broad and narrow momentum components, (b) Fit to the calculation 
of (a) with a single Gaussian, with a bin size 50 MeV/c, similar to experiment. 



